Avoiding Quantum Chaos in Quantum Computation 
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We study a one-dimensional chain of nuclear 1/2— spins in an external time-dependent magnetic 
field. This model is considered as a possible candidate for experimental realization of quantum 
computation. According to general theory of interacting particles, one of the most dangerous effects 
is quantum chaos which can destroy the stability of quantum operations. The standard viewpoint 
is that the threshold for the onset of quantum chaos due to an interaction between spins (qubits) 
strongly decreases with an increase of the number of qubits. Contrary to this opinion, we show 
that the presence of a non-homogeneous magnetic field can strongly reduce quantum chaos effects. 
We give analytical estimates which explain this effect, together with numerical data supporting our 
analysis. 
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Much attention is paid in recent years to the idea of 
quantum computation (see, for example, and ref- 

erences therein). The burst of interest to this subject is 
caused by the discovery of fast quantum algorithm for the 
factorization of integers Q demonstrating the effective- 
ness of quantum computers in comparison to the classical 
ones. Nowadays, there are different projects for the ex- 
perimental realization of quantum computers, based on 
interacting two- level systems (qubits) . One of the most 
important problems widely discussed in the literature, is 
the problem of decoherence which arises in many-qubit 
systems due to the influence of an environment [|| . How- 
ever, even in the absence of the environment, the inter- 
action between qubits may lead to the "internal decoher- 
ence" related to the onset of quantum chaos Q . 

The latter subject of quantum chaos in closed systems 
of interacting particles has been developed recently in ap- 
plication to nuclear, atomic and solid state physics (see, 
e.g., and references therein) . When the (two-body) in- 
teraction between particles exceeds the critical value, fast 
transition to chaos occurs in the Hilbert space of many- 
particle states Different aspects of this transition 
are now well understood, such as statistical description 
of eigenstates and the onset of thermalization in finite 
systems (see, e.g., [|J and references therein). 

Direct application of the quantum chaos theory to a 
simple model of quantum computer B] has shown that for 
a strong enough interaction between qubits the onset of 
quantum chaos is unavoidable. Although for L = 14— 16 
qubits the critical value J cr for quantum chaos thresh- 
old is quite large, with an increase of L it decreases as 
Jcr ~ 1/L. From the viewpoint of the standard approach 
for closed systems of interacting particles, the decrease 
of the chaos threshold with an increase of qubits looks 
generic. However, in this Letter we demonstrate that 
this conclusion is not universal and the quantum chaos 



can be avoided, for example, with a proper choice of the 
external magnetic field. 

Our consideration is based on the one-dimensional 
model of L nuclear 1/2— spins subjected to the time- 
dependent magnetic field of the following form pj , 

B(t) = [b p ^cos{v p t + (p p ),~b p L sin{v p t + ip p ),B z ], (1) 

where B z is a constant magnetic field oriented in the 
positive z-direction, 6^ , i/ p , and (p p are the amplitudes, 
frequencies and phases of a circular polarized magnetic 
field which is given by the sum of p = I, .... P rectangu- 
lar pulses of the length t p+ i — t p rotating in the (x,y) — 
plane, and providing a quantum computer protocol. The 
quantum Hamiltonian of this system has the form, 
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where the "pulse function" Q p (t) equals 1 only dur- 
ing the p-th pulse, for t p < t < t p+ i. The quan- 
tities Jk, n stand for the Ising interaction between two 
qubits , ujk are the frequencies of spin's precession in the 
B z — magnetic field, Q p is the Rabi frequency of the p- 
th pulse, If' v ' z = (l/2)a^,' v ' z , the latter being the Pauli 
matrices, and 1^ = If. ' ± il%. 

For the p-th pulse, the Hamiltonian (^|) can be repre- 
sented in the coordinate system which rotates with the 
frequency v p . Thus, for the p-th pulse, our model can be 
reduced to the stationary Hamiltonian, 
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which describes the evolution for t p < t < tp+i. 

The physically important regime of quantum compu- 
tation is a selective excitation which corresponds to the 
range of parameters: Q p <C Jk, n ^ S^k ^ ^ki where 
fefc = \u)k+i ~ w /sl IP-IT However, in this Letter, we 
consider a regime of non-selective excitation which is de- 
fined by the conditions, fl p ^> 6u>k ^> J, see details in 
[pTofl . This inequality provides the simplest way to pre- 
pare a homogeneous superposition of 2 L states needed 
for implementation of both Shor and Grover algorithms. 
In what follows we assume, for simplicity, that the am- 
plitudes of magnetic pulses, and their frequencies and 
phases are equal, tt p = Q, v v = v, and (p p = tt/2. Our 
main interest is in the nearest neighbor interaction (N- 
interaction) between qubits for two different cases, the 
dynamical one when all coupling elements are the same, 
Jk.n = J <5>i,fc+i, arid the case when all values Jk,k+i are 
random (random model). However, we will also discuss 
another case when all qubits interact to each other (A- 
interaction) with random Jk, n - In contrast to previously 
discussed model Q with homogeneous magnetic field, be- 
low we consider the magnetic field with constant gradient 
along the z-axes. Therefore, we assume that the spin fre- 
quencies Uk are slightly dependent on k in the following 
way, 8uJk <C tou- For the dynamical N— interaction and 
v = ujq, ujk = co>o + a k i a > 0); the Hamiltonian takes the 
form, 
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where 5% — ak. In the z-representation the Hamiltonian 
matrix of size N = 2 L is diagonal for = 0. For ^ 0, 
Hk,k+n — ifl/2 and H n k = H^ n . The matrix is very 
sparse, and it has quite specific structure if the basis is 
reordered according to an increase of the number s writ- 
ten in the binary representation, s = ii,Z2, ■■■■,ih (with 
i s = or 1, depending on whether a single-particle state 
of i— th qubits is the ground state or the excited state). 

The structure of eigenstates in this basis has been nu- 
merically studied in detail for L = 8, 10, 12 in dependence 
on the strength of interaction, J, and different values of 
il. For £1 ^> J the spectrum consists of L + 1 narrow 
bands with large gaps of size w fl between the bands. 
Since the most interesting energy region for preparation 
of the homogeneous wave function is in the middle of the 
energy (quasi-energy) spectrum, we consider below only 
the central band and the corresponding eigenstates for L 
even. 



In the absence of interaction between qubits, J = 0, all 
eigenstates in the given basis are fully extended with the 
value of components ip n close to 1/ \/~N. Typical structure 
of eigenstates in the central band is shown in Fig. |l] for 
different values of J. One can see that with an increase 
of J— interaction, the probabilities, w n = \ip n \ 2 , deviate 
from the unperturbed value, w n — 1/N , thus resulting 
in quantum computation errors. The data demonstrate 
the transition from regular states for a weak interaction, 
J < 0.1, to chaotic ones for J w 100. 
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FIG. 1. Typical structure 

of eigenfunctions for J = 0, 2 ■ 10~ 4 , 0.1, 1, 10, 100 denoted 
by (a,b,c,d,e,f) respectively. Eigenstates are taken from the 
central energy band for L — 12, Q — 100, a = 1. 

Global properties of eigenstates can be characterized 
by the number, N pc , of principal components, deter- 
mined through the inverse participation ratio, N pc (E) = 

E« IV'nC-E')! 4 ] ) where ip n (E) is the n-th component of 
a specific eigenf unction. For zero interaction (J = 0), 
the number N pc is equal to N, and it decreases with an 
increase of interaction, thus, giving the measure of the de- 
struction of unperturbed (J = 0) eigenstates. Note that 
for completely chaotic eigenstates with gaussian fluctua- 
tions for ip n , one has iV pc = N/3, 0. 

Numerical data for the averaged N pc in dependence on 
the interaction J define different regions, see Fig. g. The 
first region on the left of this figure corresponds to com- 
pletely "delocalized" eigenstates with ip n w l/y/~N. Here 
the energy spectrum is characterized by many close quasi- 
degenerate levels. In the second region where N pc <C N, 
all eigenstates are strongly influenced by the inter-qubit 
interaction. In this region which we denote as the re- 
gion of weak chaos, the level spacing distribution P(s) is 
quite close to the Poisson. From the data, the transition 
to the region of weak chaos occurs for J = J cr w 0.05. 
One should stress that from the practical point of view 
the weak chaos should be avoided in the quantum com- 
putation because of large deviations of eigenstates from 
the unperturbed ones, see Fig. lc-d. 

Second transition to a strong quantum chaos occurs for 
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J > 20. By the latter term we denote the situation when 
the level spacing distribution has the Wigner-Dyson form 
and fluctuations of components tp n are close to the gaus- 
sian ones with N pc m A/3. Analysis shows that this 
transition corresponds to the overlapping of the central 
energy band to the nearest ones. As one can see, strong 
quantum chaos emerges for extremely strong interaction. 
Quite unexpectedly, the results for the random model 
(with random distribution of the interaction strength in 
the interval [—J, J]) turn out to be the same. This means 
that, both for the dynamical and random A-interactions 
our system is close to an integrable one for well separated 
energy bands. 
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FIG. 2. The average number of principal components N pc 
as a function of J for the parameters of Fig. [ij; the horizon- 
tal line corresponds to N/3. The average is taken over the 
eigenfunctions from the central band only. 

The above numerical analysis for L = 12 shows that 
weak quantum chaos can significantly influence the struc- 
ture of eigenstates, but the regime of strong chaos is not 
achievable in quantum computation. However, according 
to a common believe, the threshold for both weak and 
strong chaos is expected to decrease with an increase of 
L, thus, leading to the destruction of unperturbed states 
even in the presence of a very weak interaction J. 

We show now analytically, that contrary to the stan- 
dard point of view, in this model the chaos border J cr is 
independent on the number of qubits. In order to explain 
this unexpected phenomena, it is convenient to represent 
the Hamiltonian (^) in the basis in which it is diagonal 
for non-interacting (J = 0) qubits, 



H = H a + JV . 



(5) 



Here the "effective field" Hamiltonian, Ho, is determined 
by the sum of L individual Hamiltonians H k , 
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Ho = J2 H ^ = J2 V(wfc-") 2 + n 2 I z k , (6) 



and the interaction Vq = Vdiag + Hand + Voff between 
new "quasi-particles" has the form, 



Vdi a g - -2 bkbk+\lklk+ 
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V ba 



(7) 



Voff =2E(a^W^f+i + a k +ib k I%I% +1 ) . 

k 

Three terms in the expression for Vq stand for the diag- 
onal, inside-band and between-bands interaction respec- 
tively, and 



b k = 



V - Uk 



y/(LJ k ~ Vf + ft 2 



^J(iOk - y) 2 + ft 2 
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Under the condition, \uj k — v\ <C f2, we have the follow- 
ing expression for single-particle "quasi-energies" corre- 
sponding to the Hamiltonian Hq, 



1 



1 



a 2 k 2 



(9) 



Note, that each "quasi-particle" can have 2L different 
"quasi-energies" Ck in contrast to actual particles (qubits) 
which have only two values of energy. The above expres- 
sion allows us easily to construct many-particle unper- 
turbed quasi-energies, E c — X)fc=o £fc ' ms ide the central 
band. Indeed, since we have assumed f£ 3> \u> k — v\ all 
many particles levels should have, for L even, L/2 pos- 
itive £k and L/2 negative ones. As a result, the total 
number N c i, of many-body states in the central band is 
given by N c b — L!(-|!-|!) _1 and the size can be esti- 
mate as twice the maximum energy, (AE) c b = 2Ec max ^ — 
a 2 L 2 (L- 1)/8Q. 

Now, we can estimate, for the A— interaction, the 
mean level spacing SE between those many-body states 
which are directly coupled by the interaction, see Eq. (Q) . 
The value SE can be estimated as the ratio SE = 
(AE)f/Mf where Mt w L/2 is the number of many- 
body states coupled by Viand (it is, in fact, the mean 
number per line of non-zero off-diagonal elements in the 
total Hamiltonian (0)). One should stress that the en- 
ergy range (AE) / within which these states are coupled, 
is less than the total energy width (AE) c b of the cen- 
tral band. The value of (AE) f can be estimated as the 



maximal difference between energies Ec = Y] k e k and 
Ec^ = e k of two many-body states |1) and |2) of 
-Ho, that have the coupling (l|V(, OTS d|2) different from zero. 
From the expression one finds, (AE)f = 2a 2 L/tt. As 
a result, for L 1 we have, 



SE 



M f 



2a^ 



(10) 



k=0 



k=0 



This mean spacing, SE, should be compared with the 
typical value of perturbation V — JVq S . The latter can 
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be found from Vband a s V ~ J/2. Therefore, we finally 
obtain, 



J r: 



~0~' 



(11) 



Remarkably, the threshold to chaos, J cr , does not de- 
pend on the number L of qubits. The origin of this phe- 
nomenon is that the width (AE) f and Mf both increase 
lineraly in L. This fact is entirely due to the quadratic 
growth of the energy shift in Eq. (Q) for new single- 
particles states, in dependence on k . Physically, this 
effect is related to the constant gradient of the magnetic 
field in the original Hamiltonian. 

Numerical data for the number N pc of principal com- 
ponents of eigenstates inside the central energy band in 
the new basis (where Hq is diagonal for J = 0, see Eq. 
(||)) are given in Fig. 3 in dependence on J/J cr . It is 
clearly seen that below the border, J < J cr , there is a 
scaling dependence of N pc on L and f2 which confirms our 
estimate (p"l|). On the other side, for J > J cr , the value 
of Npc saturates to its maximal value N c b/3 in correspon- 
dence with random matrix predictions . We notice that 
the transition to extended states in the Hamiltonian (||) , 
which occurs for J > J cr , corresponds to the transition 
to the weak chaos given in Fig. 2. 

In the case of random model, the data turns out to 
be similar, thus, confirming our conclusion that for the 
TV-interaction the model is close to an integrable one, in- 
dependently on whether the interaction is dynamical or 
random. One should note that, indeed, if one neglects 
the interband interaction V ff, the Hamiltonian (Q) is 
integrable, as shown in |l3|]. Thus, our analysis shows 
that the border J cr does not decrease with an increase of 
number of qubits. 
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exceed the distance ft between the nearest ones, which 
leads to strong quantum chaos. The estimate for the 
critical value, L cr , which corresponds to the band over- 
lapping in absense of the J— interaction reads as L max ~ 
(fl/a) 2 / 3 » 1. 

For the model with random interaction between all 
qubits (^-interaction) , strong quantum chaos occurs in- 
side the energy bands manifested by a Wigner-Dyson dis- 
tribution for P(s) (more details will be given elsewhere). 
In this case the estimate for the energy width (AE) j , 
within which many-body states are coupled by the two- 
body interaction, is (AE) f = L 2 /2f2, and Mf = a 2 L 2 /4. 
Therefore, we get SE = (AE) f /Mf » 2a 2 /tt which is the 
same as for the N— interaction. This is an unexpected re- 
sult since generically, the chaos border for A— interaction 
is much less than for N— interaction (J cr cx 1/L 2 , see, 
e.g. Numerical data for the A— interaction in (^) 

confirm the above result (the corresponding figure is very 
similar to Fig. 3). 

In conclusion, we have shown that, in contrast to gen- 
eral believe, the chaos border in the model of L interact- 
ing qubits does not decrease with an increase of L, in the 
presence of a strong magnetic field with constant gradient 
in the z-direction. The quantum chaos which emerges for 
a very strong interaction between qubits is irrelevant to 
quantum computation. The mechanism of strong chaos 
is shown due to the band overlapping only, and can be 
avoided even for a very large number of qubits. Moreover, 
our analysis shows that for an inhomogeneous gradient of 
magnetic field (for example, when ojj- oc fc 4 ), one can ex- 
pect that the border for chaos increases with an increase 
of L. The region of parameters for quantum computation 
with selective excitation requires additional analysis. 
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FIG. 3. The average number of principal components for 
the dynamical iV-interaction as a function of J/J cr , in the 
central band in the Ho— basis, for L = 8, 10, 12 and different 
values of Q. 

Now we have to take into account that with an in- 
crease of L, the width (AE) c b of the central band can 
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